Kondo effect of an adsorbed cobalt phthalocyanine (CoPc) molecule: 
the role of quantum interference 



o 
o 



(N 



c3 



I 

C 

O 

o 



> 

in 
in 
in 
o 

O 

c3 



c 

o 
o 



X 
S3 



G. Chiappe 1 ' 2 and E. Louis 1 

1 Departamento de Fisica Aplicada, Unidad Asociada del Consejo Superior 
de Investigaciones Cientificas and Institute) Universitario de Materiales, 
Universidad de Alicante, San Vicente del Raspeig, Alicante 03690, Spain. 
2 Departamento de Fisica J.J. Giambiagi, Facultad de Ciencias Exactas, 
Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina. 

(Dated: February 6, 2008) 

A recent experimental study showed that, distorting a CoPc molecule adsorbed on a Au(lll) 
surface, a Kondo effect is induced with a temperature higher than 200 K. We examine a model 
in which an atom with strong Coulomb repulsion (Co) is surrounded by four atoms on a square 
(molecule lobes), and two atoms above and below it representing the apex of the STM tip and an 
atom on the gold surface (all with a single, half- filled, atomic orbital). The Hamiltonian is solved 
exactly for the isolated cluster, and, after connecting the leads (STM tip and gold), the conductance 
is calculated by standard techniques. Quantum interference prevents the existence of the Kondo 
effect when the orbitals on the square do not interact (undistorted molecule) ; the Kondo resonance 
shows up after switching on that interaction. The weight of the Kondo resonance is controlled by 
the interplay of couplings to the STM tip and the gold surface, and between the molecule lobes. 

PACS numbers: 73.63.Fg, 71. 15. Mb 



Coupling of localized spins to conduction electrons may 
lead to a transport anomaly known as the Kondo effect 
P, 0. This effect, that usually shows up at low tem- 
peratures, consists of a sharp peak at the Fermi level, 
whose half-width is known as the Kondo temperature 
(Tjc), and a conductance close to one conductance quan- 
tum Qo = 2e 2 /h. The Kondo temperature in the case 
of magnetic impurities in non-magnetic metals is around 
50 K 1], whereas in artificial atoms (quantum dots) is 
just a few hundred mK 0. El In a recent experiment 
1^, |(| it has been shown that it is possible to control 
the characteristics, and even the existence, of the Kondo 
resonance by modifying the chemical surroundings of a 
magnetic atom. The experiments were carried out on 
a cobalt phthalocyanine molecule (CoPc) adsorbed on 
a Au(lll) surface. Dehydrogenation of this molecule 
(d-CoPc) by means of voltage pulses from a Scanning 
Tunneling Microscope (STM) triggered a Kondo effect 
with a rather high Kondo temperature (Tr rzs 200 K). 
This temperature is even higher than that observed for 
bare Co adsorbed on a similar surface 0, 0- Besides 
such a high Tk, one of the most remarkable results of 
[j| is the fact that the undistorted molecule does not 
show a Kondo effect, while it is readily promoted by 
distorting the molecule upon dehydrogenation. Topo- 
graphic images taken by means of the STM |5J indicated 
that the CoPc molecule has four almost non-overlapping 
lobes symmetrically placed around the Co atom. Dehy- 
drogenation distorts the molecule and forces those lobes 
to overlap. In addition it strongly decreases the distance 
from the molecule lobes to the gold surface and increases 
the Co/gold surface distance in approximately 30% 

We hereby propose a simple model that accounts for 



some of the salient features of the experiment described 
above. We take a model Hamiltonian on a small atomic 
arrangement which is solved exactly, and subsequently 
connected to semi-infinite chains used to describe the 
STM tip and the gold surface. Fig. 1 depicts this atomic 
arrangement. A central site with a single atomic or- 
bital and a strong Coulomb repulsion accounts for the 
Co atom, while the four lobes of the molecule are de- 
scribed by four atomic orbitals placed on a square whose 
center is the Co atom. Two additional orbitals located 
above and below the Co atom are included to represent 
the apex of the STM tip and an atom on the gold sur- 
face, respectively. All atomic orbitals are assumed to be 
isotropic (s-like). The Hamiltonian takes the form, 



H = 



E 



ti,3 c i<T C i<r + Un C o-\n C oi (1) 



where c i(7 creates an electron at site % with ^-component 
of the spin a —\, j, while the occupation operator ncoa 



associated to Co is, nc c 



-Coa 



CCo 



We just consider 



the half-filling case, i.e., one electron per atomic orbital. 

The parameters of this model Hamiltonian are the fol- 
lowing, tij is the hopping between atomic orbitals lo- 
cated on sites i and j (the symbol <> in Eq. (1) indicates 
that i =fi j), each orbital has an energy e,, and the local 
(Hubbard-like) Coulomb repulsion on Co is described by 
U. In particular we use the following parameters. We 
take eco — —U/2 (symmetric case) and the rest of atomic 
orbitals lying at zero energy. The hoppings incorporated 
in the model are: tc ,t (Co and the STM tip) tc ,Au (Co 
and the gold surface) tc a ,i (Co and the molecule lobes) 
(hopping between lobes). Another important param- 
eter of the model is the hopping between the lobes and 




FIG. 1: Cluster of atoms utilized to describe the CoPc 
molecule adsorbed on a Au(lll) surface. The four atoms on 
the square account for the four molecule lobes, while the atom 
at the center represents Co. The upper (lower) atom accounts 
for the apex of the tip of the STM microscope (the gold sur- 
face). To calculate the conductance semi- infinite chains were 
attached to the top atom, and to either the gold atom or the 
lobes (see text). 



the gold surface t^Au- One lead is attached to Co and 
describes the STM tip. The other lead (the Au surface) 
is attached to either Co or the lobes. 

When the cluster is connected to electrodes, the trans- 
mission across the system is given by T(E) — =j-Tr[tH] 
0, and the conductance is Q — T(Ep), where Ep 
is the Fermi level. In this expression, matrix t is 

where Ti 



r l/2 G( _ )r V2 



U(L) 



*( S U(L) 



-.(+) 



i(±) 



-ifjn-jj, ^u(L) being the self-energies of the 
upper (U) and lower (L) leads, STM tip and gold sur- 
face, respectively. Superscripts (+) and (-) stand for re- 
tarded and advanced. The Green function is written as 



G(±) = 



G 



(±) 



where G 



(±) 



is the Green function associated to the isolated cluster, 
which is obtained by exact diagonalization 0, [n| . The 
electrodes are described by means of semi-infinite chains. 
This method is exact only as far as the calculation of the 
Green function of the isolated cluster is concerned, and 
does not account for correlation effects that extend be- 
yond its bounds. The method has already been applied 
to a variety of transport problems in- 
cluding transport through hydrogenated Pt nanocontacts 

In carrying out calculations we have taken the hopping 
within the semi-infinite chains (leads) to be 1 eV, which 
is sli ghtl y larger than the hopping between s-orbitals in 
gold |13 • I n addition we take U=8 eV, not far from the 
value 6-7 eV recently estimated for the CuPc molecule 
[l7j . The rest of the model parameters have been var- 
ied aiming to identify their role in the behavior of this 
system. All calculations were done at zero temperature. 



0.5 



1 1 1 


I 1 


Ml 1 


I 1 1- 


- A , 




B 




ii 








ii 

i 1 
1 1 
i i 






- 


) 


L 




t 
t 
\ 
\ 

\ 

: 


1,1,1,1,1 


i , i , i , i , i * 



-0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 
E (eV) E(eV) 





-0.4-0.2 0.2 0.4 

E(eV) 



0.5 1 1.5 2 
t(U)(eV) 



FIG. 2: Transmission T (in units of the conductance quan- 
tum) versus the energy E (in eV) referred to the Fermi energy, 
and phase difference $ between the direct path from the gold 
surface to the STM tip going through the Co atom, and the 
path that passes through the molecule lobes (see text). A: 
ti,Au = 0, tco.t = too, a«=0.25 eV and t Co ,i=1.25 eV; tu = 
(green), and ti,i= 2 eV (red). B: tco.Au = 0, tc o ,t=0. 25 eV, 
ti t Au= 1 eV, and, i C o,i=l-25 eV; t u = (green), and t h i = 
2 eV (red). The standard Kondo resonance, obtained with 
tco,l= and tc ,t = tco, Au~0.25 eV (black), is plotted in 
both (A and B). C: phase difference calculated for the param- 
eters used in A. D: spin-spin correlation for Co/lobes (contin- 
uous line) and Co/STM tip (broken line) and total charge on 
the lobes plus Co (continuous blue line) versus ti,i, for the 
parameters used in A. 



The results depicted in Fig. 2 are addressed to identify 
the origin of the emergence of the Kondo resonance upon 
distortion of the molecule 0, 0, ■ The Figure shows 
results for either the Co atom (Fig. 2A) or the lobes 
(Fig. 2B) connected to the gold surface. In addition we 
take either non-interacting lobes (green curves) or a fi- 
nite hopping between lobes (red curves). For comparison 
we also show the results for the standard Kondo effect 
in which the Co is decoupled from the molecule lobes 
(black curves). The most appealing result is that when 
the molecule lobes do not interact no Kondo effect shows 
up [2l|. The origin of this remarkable result is likely re- 
lated to quantum interference |1 HI HI El El El l2o| . 
To illustrate this assessment we have calculated the phase 
difference between the direct path from the gold surface 
to the STM tip going through the Co atom, and the path 
that passes through the molecule lobes. This phase dif- 



ference can be derived from the following element of the 
Green function, 

G {+ \Au,t) = g(+\Au,t) + 

4g i+) (Au,Co)Y;(Co,l)G {+) (l,t) (2) 

where T,(Co,l) is a many-body self-energy that accounts 
for the lobes/Co coupling, and lower case "</' are the 
Green functions in the case that lobes and Co are de- 
coupled, (note that E(Co,/) oc t(Co,l)) 0. The results 
shown in Fig. 2C are just the phase difference between 
the two terms in the r.h.s of Eq. (2). When there is no 
hopping between lobes, the phase difference is n indicat- 
ing that the two terms may totally cancel each other as 
actually occurs j^- An alternative way to look at this 
issue is to calculate the local density of states (LDOS) on 
the Co atom that is obtained from the diagonal element 
of the Green function, 

G (+) (Co,Co) = {uj-ei)g(+\Co,Co)x 

(uj - ei )£ 2 (Co,% (+) (C ,C )]~(3) 

where uj is the energy referred to the Fermi energy. It is 
readily seen that when the lobe orbitals lie at the Fermi 
energy (e; = 0) the Green function G^ + ' {Co, Co) van- 
ishes at that energy and, thus, the LDOS at the Co atom 
[2sl Eflj | . A null density of states at the Fermi energy on 
the strongly correlated Co atom implies that no Kondo 
resonance should show up, in accordance with the phase 
analysis. Full cancellation of the two terms in the r.h.s. 
Eq. (2) is removed when the lobes orbitals do not lie 
exactly at Ep, a result that can be derived from Eq. (3). 
When coupling between lobes is switched on, the phase 
difference is no longer tt (Fig. 2C) and the Kondo reso- 
nance shows up (Fig. 2A). This results from the fact that 
switching on that coupling opens new paths for the elec- 
trons to go from the lobes to Co that contribute to the 
phase difference. Besides, the peak width becomes sig- 
nificantly widened with respect to the standard Kondo 
effect (black curve in Fig. 2A) despite of the fact that 
the molecules lobes are not connected to the gold surface. 
Fig. 2D shows the spin-spin correlation for Co/lobes and 
Co/STM tip (a similar result is obtained for Co/Au). 
Remarkably, switching on the lobe/lobe coupling shifts 
the antiferromagnetic correlation from the Co/lobes to 
the Co/STM tip (in the undistorted molecule the spin 
on the Co orbital is screened by the spin on the lobes). 
In addition, six electrons are localized in the undistorted 
molecule (lobes plus Co) a number that is reduced down 
to five when the lobe/lobe coupling is switched on (see 
Fig. 2). These results are consistent with the existence 
(absence) of a magnetic moment on the distorted (undis- 
torted) molecule, as derived from the ab initio calcula- 
tions reported in We believe that the mechanism 
hereby put forward for switching on and off the Kondo 
resonance may apply to a variety of situations. 
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FIG. 3: Transmission T (in units of the conductance quan- 
tum) versus energy E (in eV) for the Co uncoupled from 
the gold surface, i.e., tco,Au = 0, and tc a ,i = 1-25 eV, and 
tj,j = 2 eV. A: tco,t = 0.25 eV (continuous lines); ti t Au = 0.25 
eV (green) ti t Au = 0.5 eV (red) and ti,Au = 1 eV (black). 
B: tco.t = 0.5 eV (broken lines); t^Au = 0.25 eV (green) 
ti,Au = 0.5 eV (red) and ti t Au = 1 eV (black). One of the 
curves in A has also been plotted in B for the sake of com- 
parison. All energies referred to the Fermi energy. C: Peak 
weight of the Kondo resonance (in Go eV, Go being the quan- 
tum of conductance) versus the molecule lobes/gold surface 
hopping, for the two values of tc ,t of A and B (continuous 
line 0.25 eV and broken line 0.5 eV). 



Cutting out the bond between Co and the gold surface, 
and switching on those from lobes to gold, does not qual- 
itatively change these results. Again, as shown in Fig. 
2B, in the absence of lobe-lobe coupling, no Kondo effect 
shows up. We note that even though the lobes/ Au sur- 
face hopping in Fig. 2B is much higher than the Co/Au 
surface hopping used to obtain the results of Fig. 2A, the 
width of the resonance is similar and the peak height con- 
siderably smaller (compare the red curves in those two 
Figures). These results suggest that providing more ways 
to hybridize the atomic orbitals on the magnetic ion to 
the continuum sates (as may be in principle occur due to 
coupling of lobes to the gold surface) may not inevitably 
be beneficial as far as the Kondo effect is concerned. 

The effects of coupling to the leads on the width of the 
Kondo resonance are illustrated in Fig. 3. All results cor- 
respond to the set of parameters we use to describe the 
distorted molecule: finite coupling between the molecule 
lobes and no hopping between the Co atomic orbital and 
the gold surface. In addition we note that, in view of the 
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experimental information reported in |5(, it seems rea- 
sonable that the lobe/Au hopping be greater than that 
related to the Co/STM tip coupling. Both Fig. 3 A and 
Fig. 3B, show the results for fixed coupling of Co to the 
STM tip and a variable hopping between the lobes and 
the gold surface. For the lower value of the Co/STM tip 
hopping (Fig. 3A), it is noted that while the the width 
of the Kondo resonance follows the expected qualitative 
trend, it depends only weakly on coupling. The effect of 
quantum interference is demonstrated by the presence of 
a non-unitary Kondo effect characterized by a conduc- 
tance smaller than 1 (see also Fig. 2). A far more impor- 
tant effect on the width of the Kondo peak is obtained 
when the coupling between the Co atomic orbital and 
the orbital at the apex of the STM tip is increased (see 
Fig. 3B). Now the same increase in the hopping param- 
eter produces a dramatic broadening of the Kondo reso- 
nance. Fig. 3C illustrates how the weight of the Kondo 
resonance evolves with the lobe/Au coupling, for the two 
values of the Co/STM tip hopping t(Co,t) of Figs. 3A 
and 3B. For the smaller value of t(Co, t) (=0.25 eV), the 
weight increases only slightly with the lobe/Au coupling, 
saturating around 0.07 Go eV for t(l,Au) » 0.7 eV. In- 
stead, for t(Co, t) — 0.5 eV the peak weight increases 
steeply, saturating for t(l, Au) > 1 eV at around 0.5 Go 
eV. These results indicate that both couplings are equally 
important, and that, in order to increase the weight of 
the Kondo resonance, the STM tip has to get as closer 
as possible to the Co atom. We note that the results 
of Fig. 3 that more closely resemble the experimental 
data [5j, as far as the peak width is concerned, are those 
corresponding to the lower value of the t(Co,t) hopping 
(compare Fig. 3A of this work and Fig. 2 A of Ref. |5j). 

An interesting feature of our results is that upon 
switching on the hopping between the molecule lobes, 
the electron-hole symmetry is broken and the Kondo res- 
onance is no longer peaked at the Fermi level. The peak 
is displaced either upwards or downwards depending on 
the sign of that hopping. As we have assumed that all 
orbitals are s-like, all hoppings are positive. Interestingly 
enough this shifts the Kondo peak below the Fermi en- 
ergy, in agreement with the experiments 0. We finally 
note that if the calculations of Fig. 3 A are done for one 
lobe decoupled from the other three (a situation that may 
represent the case of Fig. IE of Ref. a transmission 
smaller than 0.1 Go is obtained. It would be interesting 
to check experimentally this prediction. 
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